Stability of thin liquid films and sessile droplets under confinement 
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The stability of nonvolatile thin liquid films and of sessile droplets is strongly affected by finite size 
effects. We analyze their stability within the framework of density functional theory using the sharp 
kink approximation, i.e., on the basis of an effective interface Hamiltonian. We show that finite 
size effects suppress spinodal dewetting of films because it is driven by a long- wavelength instability. 
Therefore nonvolatile films are stable if the substrate area is too small. Similarly, nonvolatile droplets 
connected to a wetting film become unstable if the substrate area is too large. This instability of a 
nonvolatile sessile droplet turns out to be equivalent to the instability of a volatile drop which can 
atttain chemical equilibrium with its vapor. 

PACS numbers: 68. 08. Be wetting, 

68.43.-h ehemisorption/physisorption: adsorbates on surfaees, 

68.03.Cd surface tension and related phenomena, 

82.60.Nh thermodynamics of nucleation 
Keywords: wetting, nanofluidies 
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I. INTRODUCTION 

Dewetting of fluid films and the ensuing formation of 
sessile droplets are both part of everyday experience. 
Moreover these mechanisms are important for the func- 
tioning of biological systems as well as for numerous 
technological processes. Dewetting on homogeneous sub- 
strates and the subsequent formation of droplets have 
been studied both experimentally [H-Q and theoretically 
p^l - [T7| in great detail. Both mechanisms can be under- 
stood quantitatively within the well established theory 
of wetting phenomena [l8l - [2l| . More recently, wetting 
and dewetting on structured surfaces receives increasing 
attention, in particular with a view on controlling the 
dewetting process on patterned surfaces [23 - |29| as well 
as in the context of microfluidics [30l - [35| . Chemical pat- 
terns consisting of lyophilic and lyophobic patches as well 
as topographic patterns such as pits and grooves effec- 
tively lead to a lateral confinement of wetting films and 
droplets. 

It is well known that confinement modifies the struc- 
tural and thermodynamic properties of condensed mat- 
ter. In small scale systems these finite size effects can 
either stabilize or destabilize certain structures. For ex- 
ample, systems exhibiting a long-wavelength instability 
are characterized by a critical wavelength such that fluc- 
tuations with larger wavelengths grow exponentially in 
time. This type of instability is suppressed in systems 
smaller than this critical wavelength. On the other hand, 
certain structures can only exist if they are larger than 
a certain critical size, such as droplets which, at least 
within classical nucleation theory, have to be larger than 
the critical nucleus. This means that certain structures 



* rauscher@is.mpg. del 
t |dietrich@is. mpg.de] 



are suppressed by finite size effects, or, to put it differ- 
ently, the availability of large space can stabilize them. 

Spinodally unstable flat films show a long-wavelength 
instability such that the dependence of their stability on 
the substrate size is obvious. Droplets of nonvolatile flu- 
ids, however, are usually considered to be stable. But 
they are in chemical equilibrium with an adsorbate or a 
wetting film connected to them [l8l - l2l| which, on very 
large substrates, acts like a reservoir: a spherical droplet 
of 100 nm radius has the same volume as an adsorbate 
layer with an effective thickness of 1 A on a substrate of 
6.5 X 6.5 /im^. Therefore, an isolated droplet of a non- 
volatile fluid placed on a macroscopically large substrate 
is expected to be unstable with respect to the formation 
of a film. If a droplet is volatile, i.e., in chemical equilib- 
rium with its vapor, it is expected to be unstable, too, 
but with respect to evaporation or condensation and the 
formation of an equilibrium wetting layer. 

The Laplace pressure in droplets decreases upon in- 
creasing their diameter while the pressure in wetting films 
is determined by the disjoining pressure. In the case of 
stable films it increases with their thickness. In a sta- 
tionary situation the pressure in the droplet is balanced 
by the pressure in the connected film. Moving a small 
amount of fluid from a droplet into its attached film in- 
creases the pressure in the drop and, as the thickness of 
the fllm increases, also the pressure in the film. However, 
due to volume conservation, the larger the substrate the 
smaller is the increase of the ensuing film thickness and 
therefore the smaller is the increase of pressure in the 
fllm. This implies, that beyond a certain substrate size 
the pressure increase in the drop is larger than the pres- 
sure increase in the film and the drop will dissipate into 
the large fllm [s^ . 

On the other hand, a substrate of limited size can only 
support droplets with a base radius smaller than half the 
substrate diameter. This means, that one can expect 
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that there is a window of droplet sizes for stable droplets 
as shown for two-dimensional droplets with small slopes 
(i.e., liquid ridges with a small contact angle) in Ref. j36| . 
Since droplet volumes scale with the third power of the 
droplet radius while the volume of the wetting or adsor- 
bate film scales with the second power of the substrate di- 
ameter, the influence of the wetting or adsorbatc film on 
droplet stability is most important on the nanoscalc be- 
cause in this case the volumes of the liquid in the droplet 
and in the film are comparable. In addition, due to the 
non- vanishing width of the three-phase-contact line there 
is a minimal size for well defined droplets [lol . [STj , which 
gives rise to an additional contribution to the finite size 
effects. 

In this spirit, here we study the infiuence of substrate 
size on the stability of fiat films and of droplets using the 
framework of density functional theory within the sharp 
kink approximation, i.e., by minimizing the correspond- 
ing effective interface Hamiltonian 38 1 in the presence of 
an effective interface potential [s^, 40 1 . 



II. EFFECTIVE INTERFACE HAMILTONIAN 

Within the capillary model for nonvolatile fiuids [4ll - 
l43j interfaces and contact lines are geometrical objects 
of zero volume and area, respectively, and the free en- 
ergy of a fluid in contact with a substrate is given by 
bulk, interface, and line contributions which are propor- 
tional to the volume, interface areas, and contact line 
lengths, respectively. Within this macroscopic model, fi- 
nite size effects occur only if the three-phase-contact line 
of a droplet reaches the lateral boundary of the substrate. 
Wetting transitions and the dependence on temperature 
and pressure of the thickness of wetting layers cannot be 
described within this macroscopic model. 

For this reason, in order to access mesoscopic scales, 
we resort to the effective interface model as the sim- 
plest non-trivial model to describe a fiuid in contact 
with a substrate. It can be derived from a classical den- 
sity functional theory using the so-called sharp-kink ap- 
proximation [3^|4^. As in the capillary model, also in 
this approach interfaces are only two-dimensional mani- 
folds but contact lines, such as the three-phase-contact 
line between fluid, vapor, and substrate have a nonzero 
width as a result of explicitly taking into account the fl- 
nite range of intermolecular interactions (for reviews see 
Refs. [lOiim). Accordingly, within this model line ten- 
sions emerge and are not input parameters psj - lsij . 

The effective, local interface Hamiltonian H for a liquid 
film in Monge parameterization z = h{x, y) on a homo- 
geneous substrate with the substrate-liquid interface A 
located in the xy-plane reads 



dx dy 



+ {d^hf + {dyhf 



with the liquid-gas interface tension a. 4>{z) is the ef- 
fective interface potential [H, IH, H^l and it describes 
the effective interaction between the liquid-vapor inter- 
face and the liquid-substrate interface. The last term 
(5/i = Ap A/i is the product of the undersaturation 
A/i = yUcooxistoncc(2^) — /i at temperature T and the num- 
ber density difference Ap = piiquid — Pvapor between the 
coexisting phases, and thus it measures the thermody- 
namic distance from the bulk two-phase coexistence line. 
Within mean field theory the equilibrium configuration 
of the liquid- vapor interface minimizes 'H[h]. 

The Monge parameterization is restricted to single val- 
ued interface configurations z = h{x,y) so that droplets 
with contact angles larger than 90° cannot be described 
this way. Therefore we rewrite Eq. (jlj in a parameter 
free form also used in the finite clement code employed 
below. For arbitrary paramcterizations of the liquid-gas 
interface r(u, v) = {x{u, w), y(u, u), z{u, v)) the area of the 
surface element is dA = §^ x ^ dudv = Ghdudv with 
the interface normal vector n pointing into the gas phase 
and G = I ^ X ^ I . In the Monge parameterization this 

reduces to G = i/l -I- | V/i(r)p, i.e., the first term in the 
square brackets in Eq. ^ . The effective interface Hamil- 
tonian can be written in terms of an integral over the 
liquid-vapor interface S 



V-ir] ^ / dA- {an{u,v) + [4>{z{u,v)) + 6fiz{u,v)] e^} , 
Js 

_ (2) 

with §2 as the normal vector of the substrate-liquid in- 
terface pointing into the liquid phase, i.e., in z-direction. 

The existence of a classical density functional has been 
proven for grand canonical ensembles [s^l- Nonetheless 
the functional in Eq. ([T]) has been used successfully to de- 
scribe also equilibrium shapes of nonvolatile fluids (i.e., 
in the canonical ensemble) by fixing the liquid volume V 
via a Lagrange multiplier p. In this case, Sfi is not an in- 
dependent parameter. It turns out, that upon adding the 
constant term Sfi V (which is independent of the droplet 
shape) to the functional in Eq. ([l}, Sfi and p multiply 
the same terms such that Sp, can be absorbed into p. 
It will turn out later (see Eq. (O) that p is the pres- 
sure difference between the liquid and the vapor, and for 
droplets one has p > 0, given the choice of sign for the 
Lagrange multiplier contribution as in Eq. ([3]). Since in 
a nonvolatile system the liquid and the vapor are not in 
thermodynamic equilibrium, the pressures do not have to 
be equal. This leads to a variation principle for the equi- 
librium shape of the liquid- vapor interface of nonvolatile 
fluids. The equilibrium shape r^'"^^ (u, v) minimizes the 
functional ( dA ■ e~ = A) 



J^[r{u,v)]^ J dA • |fTn(u, w)- 



(p{h) +6pLh 



, (1) 



{z{u,v)) -p [ z{u,v) - ^ 



(3) 



In the case of the laterally homogeneous substrates 
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FIG. 1. The effective interface potential (f>{z) (full red line) 
according to Eq. (|4]) and the corresponding disjoining pressure 
Il{z) = —<j}'{z) (dashed blue line) in units of (po and (po/ho, 
respectively. The positions of the minimum of (/>(z) at z — 
ho and of its inflection point at z = hi = v^/io ~ 1.2 /lo 
are indicated by vertical dotted lines. Also shown is (f)" [z) 
(dash-dotted green line) which appears in the second variation 
operator Oh in Eq. (|10|) and which determines the stability of 
flat film solutions (see Eq. (|17[) ). 

considered in this paper, the effective interface poten- 
tial (j){z) does not explicitly depend on the lateral coor- 
dinates {x, y). However, due to the formation of droplets 
one can still find non-trivial solutions to the minimiza- 
tion problem in Eq. ([3]). The structure of (j){z) depends 
on the types of intermolecular interactions involved. The 
simplest effective interface potential for long-ranged dis- 
persion forces (described by Lennard- Jones type interac- 
tions) and at temperatures below the wetting tempera- 
ture has the form 

The potential has a minimum of depth — a-t z = ft-o 
and an inflection point at z = hi = \/3ho ~ 1.2 ho (see 
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of the unperturbed surface and (/>'(/i) denoting the deriva- 
tive of the effective interface potential with respect to the 
local film thickness. The Euler-Lagrangc equation corre- 
sponding to the vanishing of J*-^^ is 



2aHh + n(h)+p = (7) 



Fig. [J). The potential is negative for z > •y/1/2 /iq and 
approaches zero from below for z oo. The shape of 
4>{z) corresponds to that of a continuous wetting transi- 
tion 

A. Minimizing the free energy functional 

Within mean-field theory the minimum of the effective 
interface functional J- containing the volume constraint 
(Eq. ([3])) renders the interfacial free energy for the cor- 
responding stable equilibrium configuration. 

The functional in Eq. ([3]) can be minimized numeri- 
cally by means of an adaptive finite element algorithm 
implemented by the software Surface Evolver [sst . 
Therein, the liquid-vapor interface is represented by a 
mesh of oriented triangles and, by means of a gradient 
projection method, itcratively evolves towards the config- 
uration of minimal J-' (for an example see Fig. [5]) below. 

B. Variations of the effective interface Hamiltonian 

Within the framework of variational calculus, a stable 
equilibrium profile corresponds to a vanishing first varia- 
tion and a negative second variation of the functional T. 
In order to calculate them we return to the Monge param- 
eterization and introduce the perturbed interface config- 
uration z = h{x, y) with h{x, y) = h{x, y) + e ^{x, y) and 
p = p + eip, where < e <C 1 is a small dimensionless pa- 
rameter. It is straightforward to show that the first vari- 
ation <5(i) J- of J-( [h] ,p) = T{ [h] ,p) + e S^^^T + J'^) J" + 
0{t^) with respect to the interface configuration is given 
by [h = h{x,y)) 

(5^1) J" = / dxdyif [-2(7 H,, + (j)'{h)- p] 

J A 

+ ij J^dxdy (f^-^^' (5) 

with the mean curvature 



(6) 

I 

together with 

V = [ dxdyh, (8) 

J A 

where Tl{h) = —(j)'{h) is the disjoining pressure, which 
describes the effective interaction between the substrate 
surface and the film surface, and 2 a is the Laplace 
pressure, which follows from the interface tension of the 



{dlh) [1 + (dy)^] - 2 {d^h) (dyh) (d^dyh) + (d^yh) [l + [d^h^] 
tlh = ^ 

2^iTW¥T(W 
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fluid surface. For equilibrium interface configurations the 
sum of the disjoining pressure and of the Laplace pressure 
is constant. The variation with respect to the Lagrange 
multiplier p leads to the volume constraint (see Eq. ([5])). 

The second variation d^'^^F of F with respect to the 
film height can be written form quadratic in the 
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and with the second derivative 4>"{h) of the effective in- 
terface potential. For the model potential given in Eq. ^ 
(j)"{h) is shown in Fig. [T] It is positive for small h and 
negative for large h. The second variation with respect to 
the Lagrange multiplier is identical to zero. The mixed 
variation with respect to p and h leads to the second term 
in Eq. ([9]) which due to "0 = const, vanishes for perturba- 
tions "^{x, y) which conserve the volume. The stability 
of a solution of the Euler-Lagrange equation ([7]) is deter- 
mined by the spectrum of eigenvalues of Oh- A solution 
is linearly stable if all eigenvalues are positive. 

III. THIN FILMS AND NANO-DROPLETS 

On a chemically homogeneous substrate with an area 
A there exist two distinct classes of solutions of the Euler- 
Lagrange equation ([7]). One consists of /lat films with 

h^'''i^x,y) ^hf = V/A. (11) 

The other class consists of nontrivial droplet solutions 
with one or many droplets smoothly connected to a wet- 
ting film. Here we focus on solutions with a single droplet 
because in general two or more droplets connected via a 
wetting film are unstable with respect to coarsening. In 
the following we discuss the stability of flat films and 
such droplets as a function of the substrate area A, of 
the excess liquid volume 

V,^ = V-Aho^{hf~ho)A, (12) 

and of material properties encoded in (l){h)/a. 

A. Flat films 

For flat fllms with homogeneous thickness h / the Euler- 
Lagrange equation ([7]) reduces to 

p + n(/i/) = 0. (13) 

This means that for any size of the substrate area a homo- 
geneous flat film obeying Eq. (|13p is obviously a solution 



perturbation 4': 

5^'^^F = j dxdy {^Oh-^ + (9) 
with the self-adjoined operator 



(10) 
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of the Euler-Lagrange equation. It represents either a lo- 
cal maximum, a local minimum, or a saddle point of the 
free energy functional in Eq. ([3]). The curvature of the 
interface is zero and thus the liquid gas interface tension 
drops out. If /i/ minimizes the effective interface poten- 
tial one has p = (assuming that 4){h) is differentiable) . 
For a flat interface the operator Oh in Eq. pUj) . which 
determines the linear stability of the flat fllm solution, 
reduces to 

dh^-a{dl+c^) + ^"{hj). (14) 

The corresponding eigenvalue problem has the form of 
a stationary single particle Schrodinger equation with a 
potential which is constant across the domain of the sub- 
strate. In Fig. [T]0"(z) is shown for the model potential 
from Eq. ([4]). The inverse surface tension plays the role 
of the mass. 

The eigenvalue spectrum of this operator depends on 
the shape of the domain and on the boundary conditions 
at its borders. Boundary conditions corresponding to ac- 
tual substrates of finite size are rarely compatible with 
a fiat film solution because usually there is a bending of 
the interface at the edge of the domain. For example, at 
the edge of a lyophilic patch on a lyophobic substrate the 
film thickness will go to zero (or at least to a microscop- 
ically small value) and at the brim of a flat piece of sub- 
strate the fluid fllm either continues onto the side walls 
or ends with thickness zero. The two simplest types of 
mathematical boundary conditions, which allow for flat 
film solutions, are either periodic boundary conditions or 
a Neumann type boundary condition which corresponds 
to zero slope of the fllm surface at the domain bound- 
aries. The latter would correspond to upright side walls 
with an equilibrium wetting angle of 90° at a pit-shaped 
substrate. However, even for such a setup, the interplay 
of the long-ranged forces from the substrate and from 
the side wall would lead to a bending of the film surface 

[Hill. 

For a square substrate with edge length L = \fA and 
with Neumann boundary conditions the eigenvalue prob- 
lem corresponding to Oh can be factorized by separating 
the variables and the eigenfunctions are given by plane 
waves. The degeneracy of the eigenfunctions character- 



Oh 



[i+(ax/i)2-i-(a„;i) 



dx 



0"(/i), 



izcd by wave vectors of equal modulus is alleviated by the 
boundary condition. Assuming the two edges of the sub- 
strate to be aligned with the x-axis and with the y-axis, 
respectively, the eigenfunctions are given by 

^„rn{x, y) (X COS x) COS y) , (15) 

with n, m e Nq. Since *(_„)„i = *„(-m) = ^'(_„)(_„i) = 
^rnn we only consider non-negative indices. Since we con- 
sider a nonvolatile system there is volume conservation, 
i.e., J^dx dy'^nm = and therefore either n or m have 
to be positive. The corresponding eigenvalues are given 

by 

Xnm = o {^f {n-^ + m2) + (16) 
Therefore the film is linearly stable, i.e., minA„,„ > 0, 

n , m 

for 

An' ^ ' 

For substrates of infinite, i.e., macroscopic, size A this 
is the case only if (j>"{hf) > 0. For the model effective 
interface potential in Eq. the latter inequality holds 
for 

hf < h,= \^ ho ~ 1.2 ho, (18) 

Since hi > ho (see Fig. [T]) films with negative excess vol- 
umes (i.e., hf < ho (see Eq. ((T2|) )) exhibit (t>"ihf ) > so 
that, according to Eq. (fT7|) . they are linearly stable for 
any substrate size A = L^. However, even for (p"{hf) < 
flat films are linearly stable if the substrate size L is below 
the critical value Lc = 2 tt •y/(T/|0"(/i/)|. This perturba- 
tion analysis does not yield any information about the 
nonlinear stability of film solutions, i.e., whether a flat 
film has a lower free energy than a droplet. 
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FIG. 2. Interface configurations ft''^''' of nano-droplets as 
obtained by numerical minimization of the functional T in 
Eq. (O based on the effective interface potential in Eq. ([4} for 
^/(tt/iq) — 100^ and under the constraint Vcx/(^fto) ~ 0.5, 
with (l)o/<J = 0.5 {6cq — 60°, left panel (a)), and t/io/o" = 0.1 
{9cq — 26°, right panel (b)). In the projected side view (bot- 
tom row), the underlying substrate of area A is indicated in 
blue. Lengths are given in units of ho- Drop heights are mea- 
sured from above the wetting film thickness /i^'*' = 1.009 /lo 
in (a) and h^'^^ — 1.012 /lo in (b). During the iterative mini- 
mization process, the mesh size of the triangulation has been 
coupled to the evolution of the interface shape in an adaptive 
way in order to optimize the spatial resolution locally. The 
lateral boundary conditions are implemented by a constraint 
on the boundary vertices, such that their lateral coordinates 
are fixed during the minimization process while the perpen- 
dicular height coordinate can evolve freely, effectively corre- 
sponding to neutral wetting (contact angle 90°) at vertical 
side walls (not shown) or Neumann boundary conditions. 



of the attached wetting film thus forming a contact angle 
[5l| . For the systems studied here, this contact angle is 
smaller than ^^q. 

The wetting film surrounding the nano-droplet is al- 
most flat, i.e., 2a\Hh\ < |n(/i)| (see Eq. 0). According 
to this Euler-Lagrange equation ([7]), the spatially con- 
stant pressure p is approximately given by 



B. The nano-droplet configuration 

For a given area A and a certain ratio Vex/iA ho) (see 
Eq. (fT2l) ). nano-droplets with a nonzero pressure p > 
minimize the free energy T in Eq. ([3]). This is due to 
the interplay of the surface free energy densities and the 
effective interface potential, in combination with the non- 
volatility of the liquid and the finite area A of the solid- 
liquid interface. Since the difference between the liquid- 
substrate and the gas-substrate surface tensions is given 
hy a + 4>{ho) Young's law [sl] reads [l9| : 

cos6loq = l-— . (19) 

^eq denotes the equilibrium contact angle of a macro- 
scopic drop. The influence of the ratio ^o/c on the shape 
of a nano-droplet is shown in Fig. [2] A suitable deflnition 
of the contact angle of a nano-droplet is to determine the 
curvature of its surface at the apex, to inscribe the corre- 
sponding cap of a sphere which intersects the asymptote 



p^-Ii{h^^°^~>), (20) 



and thus h^ > ho implies p > (sec Fig. [T]). 

The height hw'^'' of the wetting film, the pressure p, the 
disjoining pressure Il{hw'^'^) of the wetting film, and the 
ratio between the drop free energy J-drop and the free en- 
ergy .Ffiim of a fiat film with the same excess volume are 
shown in Table [1] for several values of Vox • For decreas- 
ing values of Vex/ {A ho) with constant A, p increases. 
This is mainly due to the increasing curvature of the 
liquid-vapor interface. For the same reason the pressure 
in macroscopic drops also increases with decreasing vol- 
ume. While the free energy Tdmp of large drops turns 
out to be smaller than the free energy J-fiim of a flat fllm 
with the same excess volume, the situation is reversed for 
Vex/ {A ho) < 0.06 (the critical excess volume lies between 
0.05 ^ft-o and 0.06 A ho). This means that nano-droplets 
below a certain size become mctastable or unstable. 



6 



Vox/(A/lo) 


h'^^^^ho - 1 


pho/a 




drop/ -^film 


0.05 


0.0263 


0.1777 


0.1779 


1.0003 


0.06 


0.0219 


0.1526 


0.1522 


0.9975 


0.10 


0.0166 


0.1191 


0.1192 


0.9782 


0.20 


0.0125 


0.0922 


0.0922 


0.9130 


0.50 


0.0091 


0.0674 


0.0684 


0.7730 



TABLE I. Characteristics of nano-droplets obtained via nu- 
merical minimization of J" in Eq. (|3]) based on Eq. Q for 
A/{nhl) = 100^ and (j)o/a = 0.5; /iL°''' is the height of the 
film at the edge of the numerical domain. The pressure p 
as the value of the Lagrange multiplier for fixing the volume 
and the disjoining pressure n(/ij^''') of the wetting film sur- 
rounding the droplet as calculated from the numerically deter- 
mined hw^^ are balanced according to Eq. (|20|l . Accordingly, 
the differences between the third and fourth column indicate 
the level of numerical accuracy. J^drop/-^flim is the ratio of the 
(mean-field) surface free energy of a nano-droplet and the free 
energy of a flat film with a height hf — Vex/ A + ho- 









,(eq) 
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FIG. 3. Vertical cut through the apex of a fully numerically 
obtained interface profile (red squares) and the corre- 

sponding approximate profile /i''^''' (full blue line) consisting 
of a spherical cap resting on a flat fllm. The proflles corre- 
spond to A/{TThl) = 100^ Vcx/{Aho) = 0.5, and (^o/cr = 0.5, 
which are the parameters corresponding to Fig. [21 a) and to 
the bottom line in Table H] Although the domain for the nu- 
merical calculation is rectangular, the droplet shape is to a 
good approximation radially symmetric (r = ^ -\- y'^). The 
wetting fllm height for both profiles is h^'^^ = 1.01 ho and the 
contact angle is Qcq — 55° compared with O^q = 60° for the 
corresponding macroscopic drop. The free energies for these 
profiles agree up to the third digit. 

C. Morphological transition 

In order to analyze the morphological phase transition 
between nano-droplets and flat films as indicated by the 
numerical data discussed above, we minimize the effective 
interface Hamiltonian in Eq. ([3]) in the subspace of in- 
terface shapes hc{x,y) describing a spherical cap sitting 
on top of a flat wetting film (see Fig. [3]). For a given total 
volume of liquid, these trial profiles are parameterized by 
the contact angle a and the wetting film height hn,. The 
latter determines the fluid volume available for the drop 
connected to the film and the contact angle determines 
the drop shape. This ansatz reduces the problem of min- 



imizing T in Eq. ([3]) to a minimization problem of the 
function 

F{a,h.^)=T[h,{x,y)] (21) 

depending on the two variables a and hw with the min- 
imum at acq and hw'^\ The corresponding minimizing 

profile is denoted by hi'^'^K In contrast to the direct, full 
numerical minimization of the free energy functional in 
Eq. the function F{a,hy^) provides also a free en- 
ergy landscape in the parameter space (a, h^,). Since for 
these two-parameter trial functions the wetting film is 
perfectly fiat, the Laplace pressure vanishes and instead 
of Eq. ([201) one has 

p = -n{h(:'^y). (22) 

In the macroscopic limit, i.e., upon increasing both A 
and Vox such that 

— 77 — >■ oo with . ™ = const (23) 

hi AhQ 

one finds hw'^'' ho for the droplet solution because 
in this limit the Laplace pressure 2 a as well as the 
disjoining pressure at the cap apex vanish. The reason 
for this is that the curvature of the droplet surface goes 
to zero if the drop size diverges and that the disjoining 
pressure vanishes for large distances from the substrate 
surface. Therefore the sum of the disjoining pressure 
and of the Laplace pressure, i.e., —p, also vanishes (see 
Eq. ([7|)). The Lagrange multiplier p does not depend 
on the position along the droplet surface and, according 
to Eq. ([22l) . the disjoining pressure on the wetting film 
is also zero. Therefore a macroscopic liquid cap with 
volume Vex is formed above the level h^'^'^ = ho where 
n(/io) = 0. The numerical minimization of F(a,h-uj) 
also yields, in this limit, agq — > ^eq with 9cq given by 
Eq. p^ . Figure [3| shows the free energy landscape 
F{a,hu]) for a large drop. All points in the parame- 
ter space {a, (/iiu — ho)/{hf — ho) < 1} correspond to 
droplet solutions (see Eq. pT|) ). i.e., a ^ 0. The fine 
(h-u, — ho)/{hf — ho) = 1 corresponds to a fiat film solu- 
tion for which F{a, h^) is independent of a because the 
volume of the droplet is zero. The global minimum of the 
free energy is located at a^q ~ Ooq = 60° and h^'^'^ ~ ho- 
The contour lines of the free energy landscape close to 
the minimum in Fig. [H are almost parallel to the a axis 
and hence shape fluctuations of the liquid cap with a 
constant cap volume are more likely than volume fluctu- 
ations, i.e., fluctuations of the wetting film height h^,. As 
shown in the inset of Fig. [3| the equilibrium angle acq ap- 
proaches the macroscopic equilibrium contact angle 9cq 
from below. 

In Fig. [4] the excess volume is chosen such that hf < hi, 
i.e.; according to Eq. (jlSp the film configuration is linearly 
stable. Nonetheless, the droplet solution is the global 
minimum of F{a,hy^). However, as shown in Fig. [H^a) 
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there is a minimal droplet size Vox below which droplets 
cannot exist: reducing the droplet size the Laplace pres- 
sure in the droplet increases until it cannot be counter- 
balanced by the negative disjoining pressure —p = Yiihf^) 
(see Eq. ([22]) ) in the film (n(z) has a minimum of finite 
depth (see Fig. [T])) and the droplet drains into the film. 
In Fig. Elja) there is also a second branch of droplet so- 
lutions which arc unstable and which have a pressure p 
intermediate between the pressure of the metastable or 
stable droplets and of the flat film. For a given value 
of T4x such a droplet solution corresponds to the saddle 
point in the two-dimensional parameter space between 
the two (local) minima given by the droplet solution 
and the flat film solution. Upon reaching the macro- 
scopic limit, the unstable droplet branch asymptotically 
approaches the flat film pressure from below (Fig.[5Ja)). 

This means, that the thickness h^'^'^ of the wetting film 
surrounding the unstable droplets approaches the thick- 
ness hf = Vc^/A -I- /iQ of the flat film solution. Therefore 
the volume inside the unstable droplets (i.e., above h^) 
decreases monotonically as the macroscopic limit is ap- 
proached. Figure [5] corresponds to Fig. 12 in Ref. [S^ 
where, however, the volume rather than the pressure is 
plotted as a function of the substrate size without dis- 
cussing the stability of the solutions. We conclude that 
in this respect in essence there is no qualitative differ- 
ence between the quasi two-dimensional ridges studied 
in Ref. [s^ and the three-dimensional systems studied 
here. 

For large excess volumes with hf > hi, according to 
Eq. (|18p . the flat film solution is linearly unstable. There- 
fore it should represent a saddle point or a maximum in 
the free energy landscape. The droplet solution should 
represent the global minimum. However, as shown in 
Fig. O^b) the flat film solution for hf = 1.25 /iq > hi is 
either stable or metastable, but not unstable within this 
only two-dimensonal parameter space considered here. 
In addition there is an unphysical branch of unstable 
droplet solutions with pressures above the pressure of 
the flat film solution. The reason for this artefact is, 
that a slightly undulated film cannot be represented 
in this two-dimensional parameter space; but spinodal 
dewetting occurs via the growth of such small perturba- 
tions. According to Eq. p7)) . the critical substrate size, 
below which the instabilit y is suppressed by the finite 
size effects, is A/ {it h^) sa 6, i.e., much smaller than 

•\/^/(7r « 11.2, the smallest substrate size for which 
the present two-dimensional parameter space analysis 
predicts the existence of droplet solutions (see Fig. [5Kb)). 
In view of this inconsistency we conclude that the results 
obtained within this approximate scheme for very small 
substrate sizes are unreliable. However, the actual sta- 
bility of droplets in the macroscopic limit is correctly 
covered within this model. 

At the morphological transition a flat film and a 
droplet of equal volume have the same free energy but 
different pressure. In the theory of thermodynamic phase 
transitions, it is common to consider transitions between 




FIG. 4. The approximate interfacial free energy F{a,hw) 
close to the macroscopic limit as described by Eq. (|23|) : 
A/{-Khl) = 10*, Kx/(A/io) = 0.06, and for 4>o/a = 0.5 corre- 
sponding to Scq = 60°. Ff = A{(T + <j){hf)) is the free energy 
of the flat film solution for these parameters; hf — ho = Ve^/A. 
The global minimum is located at q ~ 6cq = 60° and /i™ ~ ho. 
The contour lines are almost parallel to the a-axis. (The 
contour lines shown range from 0.9792 to 0.98 in steps of 
0.0001 and from 0.98 to 1.0 in steps of 0.001.) The inset 
shows the equilibrium angle Qcq upon approaching the macro- 
scopic limit as described by Eq. (|23() for the same excess vol- 
ume as used in the main figure; Qcq approaches O^q from be- 
low. The error bars are due to numerical inaccuracies. With 
hf/ho = Vc^/{Aho) + 1 = 1.06 < hi/ho = 1.2 (see Eq. ^) 
the fiat film solution with {h^ — ho)/{hf — ho) = 1 is expected 
to be linearly stable (i.e., metastable). But the number of 
data points calculated here is to small in order to be able to 
detect the corresponding free energy barrier (cf. Fig. [5] for a 
smaller substrate). 



states of different volume (or density) but equal pres- 
sure (or more general, between states with equal intensive 
state variables but distinct extensive ones). These states 
can spatially coexist with each other. However, the mor- 
phological transition between a flat fllm and a droplet is 
of a different nature in the sense that the droplet solu- 
tion and the flat film solution do not coexist with each 
other in space: the system as a whole switches from one 
solution to the other. This is not to be confused with 
the coexistence between a droplet and the wetting film 
to which it is connected. While the pressure in the wet- 
ting film and the pressure in the droplet are equal, this 
droplet configuration does not represent a bona fide ther- 
modynamic phase: its pressure changes with size whereas 
from a proper thermodynamic phase one would expect 
to be able to produce systems of different size but with 
the same pressure. In fact, Eq. ([T]) has the structure of 
a Ginzburg-Landau Hamiltonian but the potential <f>(/i) 
has its second minimum at /i — > cx). In this sense the 
droplet as a whole amounts to an interfacial region. 

The free energy landscape for finite systems with vari- 
ous excess volume ratios Vcx/{A ho) arc shown in Figs. [51- 
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FIG. 5. The pressure p = -^(liif') (see Eq. ([22}) (in units 
of p/ = — n(/i/)) calculated from the approximate free en- 
ergy in Eq. as a function of the substrate size A and 
of the excess volume Vex for (a) 14x/(^/io) = 0.04 and (b) 
Vix/(j4/io) = 0.25, i.e., for a fixed homogeneous film thick- 
ness hj = 1.04 /lo < hi and ft/ = 1.25 fto > ^i, respectively, 
and 00 /o- = 0.5. Note that for a fixed ratio V^-^j (^Aho) — Ucx 
one has Vc^/hg = nVaxx'^ with x = [A/{n hg)]^^^ . Open blue 
and full green symbols indicate metastable and stable states, 
respectively, and circles with crosses indicate unstable droplet 
solutions. Stable and metastable fiat film solutions are indi- 
cated by boxes and stable and metastable droplets by circles. 
The vertical line indicates the morphological transition be- 
tween stable films and stable droplets as obtained via numer- 
ical comparison of the corresponding two free energies. (A 
Maxwell construction for determining this transition point is 
discussed in, cf., Fig. Ilin A ^ ca corresponds to the macro- 
scopic limit, (a) For hf < hi, as in the present case corre- 
sponding to Vc^/{Aho) = 0.04, the fiat film solution is stable 
or metastable for all A. Droplets (lowest branch) occur for 
\/A/{tv hi) > 140 and they are stable for ^^/(Tr/ig) > 180. 
(b) Within the present free energy approximation, the fiat 
film solution is stable or metastable (although it should be 
unstable according to Eq. HHJ) and there is an unphysical 
unstable branch of droplet solutions (top branch). Droplets 
are stable or metastable for for all substrate sizes A. 



m For the large excess volume in Fig. |6l the droplet con- 
figuration with Q!oq w 57° and hw'^^ — /lo ~ 0.2 {hf ~ ha) 
is the global minimum. The flat film solution with 
hf = lAhg < hi is linearly stable as expected for the cho- 
sen effective interface potential (see Eq. (|T8)) ). Upon de- 
creasing the excess volume the free energy of the droplet 
solution increases and the minimum becomes shallower 
(see Fig. [7]). At a certain excess volume, the flat film 
solution becomes the stable solution and the droplet so- 
lution becomes metastable. Reducing the excess volume 
even further, the free energy minimum corresponding to 
a droplet solution becomes more and more shallow un- 
til it finally merges with the corresponding saddle point 
(see Fig. [8]) , and vanishes completely. This leaves the film 
solution as the only stable solution. 

This morphological transition is visualized even bet- 
ter by forming vertical cuts of the free energy landscape 
at fixed h^, i.e., parallel to the a-axis and by seeking 
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FIG. 6. The approximate interfacial free energy F{a,hu,) 
for A/{Tihl) = 100^ Kx/(A/io) = 0.10, and cjio/o = 0.5. 
Ff = A ((7 + <j){hf)) is the free energy of the fiat film so- 
lution for these parameters; hf — ho = Vc^/A. The global 
minimum representing a stable nano-droplet is located at 
(ft-m''' — ho)/{hf — ho) ~ 0.2 and Ofeq ~ 57°, i.e., close to 
but smaller than the macroscopic equilibrium contact angle 
6cq ~ 60° for this system. The fiat film solution /i^''' = hf 
(so that {hw'^'' — ho)/ {hf — ho) = 1) is metastable; for this so- 
lution there is no dependence on a. Contour lines are shown 
in the range 0.9802 to 0.9809 in steps of 0.0001 and from 0.981 
to 1.005 in steps of 0.001. 
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FIG. 7. The approximate free energy F{a,hw) as defined in 
Eq. ^ for A/{7vhl) = 100^ Vex/(^/io) = 0.06, and (j^o/a = 
0.5 (i.e., for the same parameters as in Fig.|6]but for a smaller 
value of Vox). Ff = A{a + (j){hf)) is the free energy of the 
flat film solution for these parameters; hf — ho = Vcyi/A. The 
contact angle corresponding to the global minimum is Qcq ~ 
55° at (ftm'''' — ho)/{hf — ho) « 0.04, i.e., smaller than in 
Fig. H The fiat film solution h^^'^'^ = hf (so that (/iL'"'' - 
ho)/{hf — ho) — 1) is metastable and exhibits no dependence 
on a. Contour lines are shown in the range 0.99902 to 0.9991 
in steps of 0.00002 and from 0.9992 to 1.003 in steps of 0.0002. 



the minimum of the free energy F within each cut as a 
function of a. This renders i^min(/i«)) = Tohia F{a, h^). 
In Fig. [9] the corresponding the minimal free energy 
^min(^tii) is shown as a function of the wetting film thick- 
ness h^. The energy scale is normalized by the free en- 
ergy Ff{Vex) of the corresponding flat film solution (com- 
pare Figs. [4] to [8|). For very small excess volumes the free 
energy as a function of the wetting film thickness is mono- 
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FIG. 8. The approximate interfacial free energy F{a, h^) for 
A/{-Khl) = 100^ V^^/{Aho) = 0.048, and 4>o/a = 0.5 (i.e., 
for the same parameters as in Figs. [6] and [7] but for a smaller 
value of Vex). Ff = A {a + <j){hf)) is the free energy of the flat 
film solution for these parameters; hf — ho = Vc^/A. For this 
excess volume the droplet solution has disappeared and the 
flat film solution h^'^^ ~ (^o that (/il^^'^ — fto) / (hf — ho) = 1) 
is the global minimum. Contour lines are shown in the range 
0.999 to 1.006 in steps of 0.00025. 
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FIG. 9. The minimal free energy Fminihm) ~ min^ ^"(0!, h^) 
as a function of the wetting film thickness /im for A/{TvhQ) = 
lOO'^, </)o/o" = 0.5, and several values of Vex/{Aho) ranging 
from 0.048 to 0.3. Ff = A{a + 4i{hf)) is the free energy 
of the corresponding flat film solution for these parameters. 
The morphological transition between flat films and nano- 
droplets occurs between Vcx/{Aho) — 0.058 and 0.056. For 
Vex/ {A ho) > 0.2 the flat film solution appears to becomes un- 
stable as expected from Eq. (|18p (in the inset see the enlarged 
view of the region near h^ ~ hf). The symbols indicate the 
points calculated numerically. For Vex/{Aho) = 0.3 a small 
barrier cannot be ruled out on the basis of the available nu- 
merical data. 



tonically decreasing and the only minimum which occurs 
is the one corresponding to a flat film of thickness hf so 
that {hw — ho)/{hf — ho) = 1. For intermediate excess 
volumes (0.05 < Ve^/{Aha) < 0.058 in Fig. HI) there is a 
second minimum corresponding to a metastablc droplet. 
With increasing Vex this droplet minimum deepens un- 
til it is as deep as the minimum corresponding to the 
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FIG. 10. The pressure p in units of a /ho as a function of 
Vex/ {A ho) for a fixed substrate size, A/{nhl) = 100^ and 
(/!>o/cr — 0.5 as obtained from the approximate free energy 
expression F{a,hw). Global minima (full green), local min- 
ima (open blue), and saddle points or local maxima (symbols 
with red crosses) are shown. The upper branch corresponds 
to flat fllms (boxes) and the lower one to nano-droplets (cir- 
cles). Within the reduced model we have p = —Il{h'^) (see 
Eq. (|22|) '). The pressure values obtained from a numerical 
minimization of Eq. ((3]) (black diamonds; p « —Tl{h'^), see 
Eq. (|20|) ) agree well with the results obtained from the ap- 
proximate free energy. The dashed vertical line indicates the 
volume at which hf = hi. At this volume the unstable droplet 
branch merges with the flat film branch. The full vertical line 
indicates the morphological transition between the film and 
the droplet configurations. 



flat film (at Vcx/{Aho) w 0.057). This marks the point 
of the morphological transition between a flat film and 
a droplet solution. Increasing the excess volume even 
further the droplet solution becomes more stable. Ac- 
cording to the inset of Fig. IH] it seems that the flat film 
solution (i.e., {hw — /io)/(/i/ — ho) = 1) becomes unstable 
for Vox/ {A ho) > 0.2 as expected from Eq. pB]) . However, 
for this latter value a tiny free energy barrier cannot be 
ruled out on the basis of the available numerical data. 

Figure (TU] shows the pressure p as a function of the 
excess volume for a homogeneous film of thickness hf 
(upper curve) and for the droplet solution (lower curve). 
The upper branch is exact while the lower branch is cal- 
culated by minimizing the approximate expression for the 
free energy F{a,hw) defined in Eq. ((2T|) . According to 
Eq. ([m, for both branches one has p = — n(/i^^). Fig- 
ure [TO] also shows pressure values obtained by numerical 
minimization of the full functional J- (for which, accord- 
ing to Eq. p ss -n(/i^9)). The pressure in the 
flat films (upper curve) is given by p = —Il(hf) with 
hf/ho = Vox/ A + 1 (see Eq. dH])) and it has a maxi- 
mum at Vcx/{Aho) = 0.2, corresponding to hf = hi. For 
excess volumes smaller than Vex /{A ho) = 0.2 the flat 
film solution is metastable or stable. For larger excess 
volumes, the flat film solution is linearly unstable. How- 
ever, the spinodal wavelength is extremely large close to 
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the pressure maximum such, that, according to Eq. p7)) 
and for the given substrate size, the instability actu- 
ally sets in only for Vc^/{Aho) > 0.2011. In Fig. \M 
for 0.05 < Vex/ {A ho) < 0.2 there are two curves be- 
low the curve corresponding to the flat film solution; 
the upper one (red circles with crosses) corresponds to 
a saddle point in the free energy landscape and the 
lower one corresponds to a (potentially local) minimum. 
Both branches represent droplet solutions. The unsta- 
ble branch ends at Vex/{AhQ) = 0.2, i.e., at the maxi- 
mum of the pressure in the flat film solution. The three 
curves in Fig. [TU] form a hysteresis loop. The value of 
Vex/iAho), at which the transition (thin vertical line in 
Fig. [T0| between a flat film and a droplet occurs, can 
be obtained either by comparing free energies directly 
or via a Maxwell construction (see Fig. [TT]) . The latter 
can be shown by integrating p = dJ-/dV = dJ-/dVcx 
(due to Eq. dS]) and since dV/dVex = 1 due to Eq. ([T2|) ) 
along p(y): J-{V) ~ TiVo) = Jy^piV)dV. Starting the 
integration at the volume Vcq at which the free energy 
of the film (upper branch) and the stable droplet (low- 
est branch) are equal (see Fig. [TT|) one integrates up to 
V = Vi, i.e., the volume of a film of thickness hi at which 
the unstable droplet branch merges with the film branch. 
The result, i.e., the sum of area (1) and (2) in Fig. [11] is 
the diflterence of the free energies of a film with volume 
Vi and a film with volume T4q. At Vi one switches to the 
unstable droplet branch and integrates down to its end at 
Vm- The result is the difference between area (1) and the 
sum of area (3) and area (4). From there one continues 
on the metastable droplet branch up to Vcq, which adds 
area (4). As a result, the difference of the free energy of 
a flat film of volume V^q and a stable droplet of the same 
volume is the difference between area (1) and area (3). 
For the chosen model interface potential in Eq. Q the 
flat film solution becomes linearly unstable at the value 
of Vex/{Aho) (i.e., 0.2 in Fig. [TOl), where the unstable 
droplet curve merges with the fiat film curve. 

In Fig.[Tn]the excess volume is expressed in terms of the 
substrate area. In order to discuss whether the minimal 
droplet size is determined by the interface potential or 
by the substrate size, one could fix the excess volume 
Vex (as a measure for the droplet size) and the substrate 
potential and vary the substrate size A. But the excess 
volume is defined as the fiuid volume above the height 
ho (see Eq. p^ l and increasing the substrate area A 
for fixed Vex means effectively reducing the droplet size. 
The droplet volume Vd = V — Ah^ above the height 
of the wetting film h^ is a more suitable measure for the 
droplet size. For this reason in Fig.[T2]wc plot the droplet 
volume Vd as a function of the excess volume Vex for two 
substrate sizes. The data are obtained in the following 
way: for each fixed value of A and of Vox (i-e., for fixed 
total volume V = Vex + ^ ^o) the interfacial free energies 
as shown in Figs.|6H8]are calculated. The position of local 
and global minima and of saddle points (corresponding 
to stable, metastable, and unstable droplet or flat film 
solutions) are determined numerically, in particular the 
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FIG. 11. Sketch of the Maxwell construction leading to the 
position of the thin full vertical line in Fig. 1101 The color 
code corresponds to the one in Fig. 1101 green, blue, and red 
indicate stable, metastable, and unstable states, respectively. 
Vb denotes the volume of a film of thickness /lo, Vm is the 
minimal volume required to form a droplet, Veq is the volume 
at which the free energies of the flat film and of the stable 
droplet are equal. For V /"Vi the branch of metastable flat 
flims turns into a branch of unstable flat films. There also the 
branch of unstable droplets merges into the flat film branch. 



wetting film thickness /i^, from which one can determine 
the droplet volume Vd ~ Vex — Aih^j^ho)- AsinFig. [TUJ 
for large Vex there are three branches of solutions (flat 
film solutions with = 0, unstable droplet solutions, 
and metastable or stable droplet solutions). For small 
Vex there are only flat film solutions. The size V/^ = 
Vd{Vex = Kx) of the smallest metastable droplet (which 
is identical to the size of the largest unstable droplet) 
increases with the substrate area, as well as the value V^^ 
of the corresponding excess volume. However, V^^/ {A ho) 
decreases upon increasing A (compare Figs. I12f a) and 
(b)). This means, that the thickness h'j = Vex/ A of the 
fiat film solution corresponding to the minimal droplet 
also decreases upon an increase of the substrate area. 

The nonexistence of droplet solutions for too small 
values of Vex can be rationalized by considering a fur- 
ther simplified reduced expression for the free energy. 
Neglecting the influence of the disjoining pressure on 
the spherical cap the minimization problem for F{a, h^) 
yields (see Eq. ((2T|) and up to the constant substrate- 
liquid surface tension) 

F={A-r''^)[a + 4>{hn,)]+^Sd, (24) 



with r ~ a/2 /id i? — ft,^ denoting the base radius of the 
rfrop (taken aX z = hw) and Sd = 2tt Rhd denoting the 
surface area of a spherical cap of height hd and radius 
R. The volume of the spherical cap is given by Vd = 
fhj{3R- hd) and the total fluid volume by V = Vd + 
Ahuj. It is convenient to write the volume constrained 
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FIG. 12. The droplet volume Vd = V — Ah-u, as function of 
Vc^/ [Aho) (i.e., as a function of V = \4x + Aho) for (f>o/a = 
0.5 and (a) A/ {hi tt) = 100^ and (b) A/ {hi vr) = 1000^ as ob- 
tained from the approximate expression F{a, h^) for the free 
energy: for fixed V — Vox + Aho the free energy landscapes 
(see, e.g.. Figs. [BHS} have been calculated and the wetting 
film thicknesses of the droplet solutions — if they exist — 
have been determined. The upper branch corresponds to the 
stable (green) or metastable (blue) droplet solution. For (a) 
this is the lower branch in Fig. 1101 The points on the abscissa 
correspond to stable (green) or metastable (blue) fiat film 
solutions. The comparison between (a) and (b) shows that 
the minimal droplet size (for which the unstable and the 
metastable droplet branches meet) increases upon increasing 
the substrate area while the corresponding excess volume Vox 
in units of the substrate area decreases. 
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as a function of the droplet height hd rather than the 
droplet contact angle a. The minimum of F{hd, h^u) fol- 
lows from the zeroes of its first derivatives with respect 
to hd and hw ■ Using the above expressions for V and Vd 
one obtains from dF{hd, hyj)/dhd = 



R(l){hu,) + ahd = 0. 



(26) 



Using this expression together with the above expressions 
for V and Vd one obtains from dF{hd, hw)/dhw = 0, after 
reintroducing a via the geometric condition r = R sin a, 



n(/i^) = — 



2cr 



R I 



(27) 



Apart from a small correction (which is small if A is large 
compared with the base area tt R^ sit? a of the droplet) 
Eq. (P7|) tells that the disjoining pressure in the film and 
the Laplace pressure 2 a Hh = —2(j/R (see Eq. ([5])) in 
the droplet are equal (according to Eq. (O both are equal 
to p). Using the geometric relation cos a = 1 — hd/R in 
Eq. we also get 
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FIG. 13. (a) Disjoining pressure n(/i,„) with ho < h^ < hf 
(calculated for the model potential in Eq. Q) in the wetting 
film (solid red line) and the right hand side of Eq. (I27p with 
R{K„,hf,A) from Eq. (f29)l (dashed blue lines) for A/{Tr hi) = 
10^ and (f>o jo — 0.5 as a function of the wetting film thickness 
h^ for hj = 1.04/10, hf = h} ^ 1.05005 /lo, and hf = 1.06 /lo, 
which fixes V = Ahf for a given A (see Eq. (|12l) '). We also 
show the right hand side of Eq. ((23 for A/{-Khl) = 10^ and 
hf = hf = 1.02587 /lo (dash-double-dotted cyan line), as well 
as for A/{nhl) = 10'^ and hf = hj = 1.01394 /lo (dash-dotted 
green line). The thin dotted red line shows the linear fit to 
II(/Ito) at hw ~ ho. For hf — h'f{A) the curves n(/i™) and 
the one for the right hand side of Eq. H27p touch each other 
at a single point at h^ — h'^{A) indicated by a magenta cir- 
cle, (b) hy and h'^ as a function of the substrate area A in 
units of TV hi as obtained graphically from (a) (black squares 
and magenta circles, respectively) and from the analytic ap- 
proximation described in the main text (full black and dotted 
magenta line, respectively). 



In the macroscopic limit i? — > 00 in Eq. (j27p implies 
n(/i,u,) — > 0, i.e., h^ — )- ho so that (l){hw) — > 0(/).o) = — 0o 
(see Fig. [T]), and therefore a — )■ O^q (see Eq. p^ ). As a 
function of a, R, and h^, the total conserved fluid volume 



V = Ah,, 



ttR^ 



(2 + cos a) (1 — cos a)^ 



(29) 



For a given value of a Eqs. ^F7\ and ((^ provide solutions 
for hw and R only if V is sufficiently large. The thick- 
ness hw can only vary between ho (i.e., the whole excess 
volume is concentrated in the droplet) and hf = V/A 
(i.e., there is no droplet). For h^ = ho the disjoining 
pressure 11 in the film is zero while the Laplace pressure 
2aHfi ~ —2a/R in the droplet is negative. Both be- 
come more negative for increasing /i„ because the droplet 
shrinks and n'(/io) < 0. The Laplace pressure diverges to 
—00 as hyj ^ hf = V/A because the droplet volume Vd 
(and therefore the droplet radius R) vanishes in this limit 
and Hh = — 1/i?. But the disjoining pressure is bound 
from below. With arccos[l + (t){hyj)/a] (sec Eq. (^5]) ) 
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Eq. ((29)) can be solved for R yielding R{h^, V, A) or 



R{hw, hf,A) 



3A{hf-h^) 



3 1 4>{hiu) 













(30) 



due to V = Ahf. Accordingly, one can consider both 
sides of Eq. (P7| as a function of h^^ as shown in Fig.fT^la) 
where a is approximated by O^q = arccos(l — (j)Q/cr). The 
right hand side of Eq. ([27| increases (decreases in abso- 
lute value) upon increasing hf/ho = V/{Aho). For large 
A/R"^ and sa ho, the right hand side of Eq. ([27| is 
approximately given by 



2 

?! 



(3 + ^ 

(I'D \ <^ 

TT h'i hn 



(31) 



The two curves only intersect if the fluid volume (or hf = 
V/A) is sufficiently large (see the three blue dashed curves 
in Fig. ma)). 

For sufficiently large excess volumes, i.e., for suffi- 
ciently large hf there are two intersections in Fig. [T^ a). 
Because for fixed total volume V increasing (i.e., in- 
creasing the amount of liquid in the film) means decreas- 
ing the droplet volume, the intersection at the larger val- 
ues of corresponds to the unstable solution while the 
intersection at the smaller value of hw corresponds to the 
stable droplet solution. (The unstable droplet is always 
smaller than the stable one.) In the macroscopic limit 
A/hl 00 with fixed hf ^ ho + Vcx/A (sec Eq. ((23)) ') 
the stable solution moves to h^ —> ho . This means that 
the volume of the stable droplet gets very large because 
due to hw — >■ ho the whole excess volume goes into the 
droplet. In the macroscopic limit, the unstable solution 
moves to h^ — >■ hf. We can obtain the corresponding 
leading behavior by the following procedure. First we in- 
sert R{h^, V, A) as obtained from Eq. (gni) into Eq. (|?7)) 
and we replace cos a by the expression in Eq. (|28)) . After 
substituting V = Ahf we expand both sides in powers 
oihw — hf and we obtain in leading order hf — hw ^ ^ I A 
[57| . As a consequence, in the macroscopic limit the vol- 
ume Yd = A{hf — h^) of the unstable droplet should 
converge to a finite value. However, this primitive model 
only applies to large droplet volumes and therefore this 
result for unstable drops might turn out to be an artefact 
of the approximations used. 

As shown in Fig. ITST b) the critical average film thick- 
V^/A required for forming a droplet decreases 
as a function of the substrate area. For very large A 
both h'j and the corresponding wetting film thickness h^ 
corresponding to the smallest possible droplet are very 
close to ho such that in Eq. ()27|) one can expand n(/iu,) 
around ho (see the thin dotted line in Fig. [T^a)). If one 
makes the additional approximations of using a sa 9cq 
and of reducing the right hand side of Eq. p7p to the 
Laplace pressure by neglecting the term ^ A"^ in the 
denominator of the right hand side of Eq. (|?7)) . one 



ness h'^f 



can determine /i^ and /i^, analytically with the result 
h}/w - ho (X A-^l^ (see Fig. [IStb)). The drop volume 
is Vd = V — Ahu, = A{hf — h^^). Thus the volume 
of the smallest possible droplet diverges for ^ — > 00 as 
V;[ ^ {h<) ~ hl^) A (X A^^\ 



IV. SUMMARY AND CONCLUSIONS 

We have studied the stability of nonvolatile flat films 
and droplets on smooth and chemically homogeneous 
substrates with finite surface area A. The analysis is 
based on density functional theory within the so-called 
sharp kink approximation, i.e., by minimizing the effec- 
tive local interface Hamiltonian with the effective inter- 
face potential shown in Fig. [TJ 

The stability of flat films and of nano-droplets is 
strongly affected by finite size effects. We have shown 
that in these systems (i) spinodal dewetting can occur 
only if the substrate area A is large enough to support 
the shortest unstable wavelength, (ii) there is a mini- 
mal size for droplets connected to a surrounding wetting 
layer, (iii) droplets are unstable with respect to drainage 
into a connected wetting films if the substrate area is too 
large, and (iv) that fluctuations of the droplet shape un- 
der the constraint of a fixed volume are more likely than 
volume fiuctuations. 

Our findings are manifestations of the general rule 
that long-wavelength instabilities are suppressed by fi- 
nite size effects. The shortest instability wavelength 
Lc = 2 TT \/cr/\(j)"{hf)\ of spinodal dewetting depends on 
the material properties, i.e., on the surface tension a and 
on the effective interface potential </!)(z), as well as on the 
average film thickness hf = V/A, whereas V is the con- 
served total liquid volume. In particular for film thick- 
nesses close to infiection points of (j>{z) and for thick films 
this wavelength becomes very large. For differentiable 
effective interface potentials the second derivative has a 
maximum (typically at a thickness of a few ho where 
4)'{ho) ~ 0). Therefore the spinodal wavelength Lc of 
films with the corresponding thickness has a minimum. 
Experimentally spinodal wavelengths of the order of mi- 
crons have been reported . This means that spinodal 
dewetting can be suppressed by structuring the surface, 
e.g., by a periodic pattern of hydrophilic and hydropho- 
bic stripes, the latter ones with a width smaller than Lc- 
The width of the hydrophilic stripes which is necessary 
to stabilize the film has to be determined separately. 

We have calculated the shape of nano-droplets numer- 
ically as shown in Fig. [2] and we have determined the 
thickness h^ of the wetting film on which the nano- 
droplet resides (see Table |l|. Using a subset of trial func- 
tion for the droplet shape which are parameterized by 
the contact angle of the droplet and by the wetting film 
thickness h^ (see Fig. [3]) we have mapped the free energy 
landscape of the system (see Figs. S) and [BHS])- 

In contrast to macroscopic drops (see Fig.)?]), for nano- 
droplets the influence of the wetting film to which they 
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are connected cannot be neglected. If the excess volume 
V^x = V — AhQ = {hf — Hq) a is fixed, there is a minimal 
substrate size below which no droplet solutions exist (see 
Fig. [5]). Conversely, for a fixed substrate size A one can 
find droplet solutions only above a critical (excess) vol- 
ume (see Figs. [71 and fTO]). This is reminiscent of classical 
nucleation theory which also leads to the notion of a crit- 
ical nucleus size. However, in the latter case one usually 
considers unbounded systems such that one cannot ob- 
tain stable droplet solutions at all. In the present case, 
the conserved total volume of fluid is distributed between 
a finite sized wetting film and a droplet; this allows for 
stable droplet solutions. 

As illustrated in Fig. [11] the volume V (or excess vol- 
ume = V — Aho) at which the free energy of the flat 
film solution (a film of homogeneous thickness hf = V/A) 
equals the free energy of the stable droplet (indicated by 
a thin vertical line in Fig. llOp can be determined by a 
Maxwell construction. This construction is based on the 
observation that the Lagrange multiplier p (i.e., the pres- 
sure difference between the liquid and the vapor phase) is 
given by p = dT/dV, i.e., by the partial derivative with 



respect to the chosen total volume V (see Eq. ([3])). 

The size of the smallest possible droplet increases (see 
Fig. ll2p and the thickness of the wetting film surrounding 
the droplet decreases upon increasing the substrate area 
(sec Fig. [T5)) . Within a suitable approximation of the free 
energy wc have found that the volume VJ^ of the smallest 
possible droplet diverges upon increasing the substrate 
size A as V^/Hq oc {A/hD^^^. The proportionality fac- 
tor depends on the equilibrium contact angle Ocq and for 
nonzero contact angles it is of the order of unity with 
0.077 as a lower bound (realized at 9eq = 180°). For 
/iQ ~ 1 A this means that the minimal droplet volume 
on substrates of size A = Imm^, 1/im^, and (100 nm)^ 
equals that of a cube of edge length 300 nm, 10 nm, and 
1.4 nm, respectively. On the same substrate the volumes 
of the connected wetting films of thickness 1 A fit into 
cubes of an edge length of 4.6 /im, 46 nm, and 10 nm, re- 
spectively, i.e., they are much larger. Our results show 
that nonetheless the finite extent of the substrate surface 
plays a significant role for the droplet formation and the 
associated morphological phase transition. 
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